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Many experimental systems exist that possess charge-density-wave order in their ground state. 
While this order should be able to be described with models similar to those used for superconduc- 
tivity, nearly all systems have a ratio of the charge-density-wave order parameter to the transition 
temperature that is too high for conventional theories. Recent work explained how this can happen 
in harmonic systems, but when the lattice distortion gets large, anharmonic effects must play an 
increasingly important role. Here we study the gap ratio for anharmonic charge-density wave sys- 
tems to see whether the low-temperature properties possess universality as was seen previously in 
the transition temperature and to see whether the explanation for the large gap ratios survives for 
anharmonic systems as well. 
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I. INTRODUCTION 

The appearance of stripe ordering in materials re=. 
lated to the high-temperature superconducting oxidesEl 
and of charge and orbital ordering in the colossal mag- 
netoresistance materialsB, has provided a renewed inter- 
est in the physics that drives charge-density-wave or- 
der. This phenomenon has been seen in a wide variety 
of materials ranging from quasi-one-dimensional systems 
(where Peierls distortion physics is important) such as. 
the organic conductorsEl, to the di- and trichalcogenidcsia 
like 2H-TaSe2 or NbSfe, so-called A15 materialstLsuch 
as V3S1, blue bronzesB like KM0O3, cubic oxidesQ like 
Bai-^K^BiOs, and Verwey transition materialsB such as 
Fe 3 4 . 

There has been much theoretical work on this problem 
as well. The dynamical mean field theory was employed 
to solve for the charge-density-wave phase in the ordered 
states and the puzzks-of the large gap ratio was resolved 
for harmonic systemalJ. The resolution involves properly 
accounting for phonon renormalizations, for nonvanish- 
ing effects from vertex corrections, and from effects due 
to nonconstant electronic density of states (since con- 
ventional Migdal-Elisahberg approaches cannot produce 
gap ratios 2A/T C larger than about 8). The anharmonic 
problep, was also examined with dynamical mean-field 
theorylliJ and it was found that the transition temperature 
satisfied a scaling law with the wavefunction renormal- 
ization parameter for a wide range of parameter space 
(including systems that were not Fermi liquids at low 
temperatures) . 

Here we concentrate on the issue of whether or not the 
gap ratio can remain large when anharmonic effects are 
present. The scaling law for the transition temperature 
shows that anharmonicity does not have a dramatic effect 
on T c , but as the system is cooled down to T — 0, the lat- 
tice distortion becomes larger and larger generating the 
full charge-density-wave gap. Naively, we would expect 
anharmonicity to reduce the lattice distortion (relative 
to a harmonic system) because the higher powers in the 



phonon potential do not allow the phonon coordinate to 
move as far away from the origin. Hence, one expects 
that anharmonicity will generically reduce the gap ratio, 
since T c will be unaffected, but A will be reduced relative 
to the results of a harmonic system with the same value of 
the wavefunction rcnormalization parameter. We need to 
verify whether this effect occurs and determine how large 
it can be to see whether one can still explain the large 
gap ratios of charge-density-wave systems in the presence 
of anharmonic potentials. 

In Section II we introduce the model and the tech- 
niques used to solve for the gap ratio. Section III con- 
tains our results and discussion, and Section IV contains 
our conclusions. 



II. FORMALISM 

We will be investigating the-s.tatic anharmonic Holstein 
model, whose Hamiltonian isE30 

H=-^2 Ujcj a c ja + J2(9%i - m)0m + n a) 

i,j.cr i 

i i 

Standard notations are used here: c\ a (a a ) creates (de- 
stroys) an electron at lattice site i with spin er, rii a — 
cl a Ci a is the electron number operator, /2 is the chemi- 
cal potential, and xi is the phonon coordinate at lattice 
site i. We examine the static (or classical) phonon case 
here, so there is no kinetic energy of the phonon. The 
hopping of the electrons is restricted to nearest neigh- 
bors on a hybercubic lattice in d-dimensions. We take 
the limit d — > 00 and^ale tij = t* /(2%/rf) in order to 
have nontrivial resultsEij. The rescaled hopping integral 
t* determines our energy scale (t* = 1). The bare den- 
sity of states becomes a Gaussian exp(— e 2 )/y / 7r, with e 
the band energy. The local phonon is taken as a classi- 
cal variable, so it only has a spring constant k associated 
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with it. The anharmonic potential is chosen to be of 
the simplest form, a quartic term with a strength a an . 
The deformation potential (or electron-phonon interac- 
tion strength) is denoted by g and measures an energy 
per unit length. 

It is useful to shift the phonon coordinate, in order to 
see explicitly the particle-hole symmetry present in the 
harmonic model. We shift Sj — > Xi + x' with g(n) + Rx' + 
3a an x' 3 = (with (n) the average total electron filling), 
to transform the Hamiltonian into 

H = - Ujc\ a c ja + J^igxi - /x)(n iT + n a - (n» 
+ ^ K J2 X 1 +(3 an J2x 3 +a an ^2x$, (2) 



with k = R + 12a an x' 2 , Pan — 4a an x', and /i ~ u — 
gx'. It is the presence of the cubic term, when a an 

that removes the particle hole symmetrycj from the 
problem when [i = and (n) = 1, since the Hamiltonian 
is no longer unchanged under the transformation rii a — > 

1 — Tiia- and Xi — > —Xi. Note that the system still will 
possess charge-density- wave order at half filling and weak 
coupling (for all small g > 0) because the band structure 
is still nested at half filling, even though the Hamiltonian 
is not particle-hole symmetric. 

Our calculations are performed using-standard tech- 
niques of dynamical mean-field theorytEJ: we iterate a 
series of equations to self-consistency that involve (i) de- 
termining the local Green's function from the self energy 
(by integrating over the nonintcracting density of states), 
(ii) extracting the effective medium (by removing the self 
energy from the local Green's function), (iii) calculat- 
ing the probability distribution of the phonon coordinate 
(by solving the atomic path integral in a time-dependent 
field), (iv) extracting the electronic self energy (after de- 
termining the local Green's function from an integral over 
the phonon coordinate distribution). This procedure is 
standard (even for the case of the ordered phase), and 
we do not describe any further details here. Our calcu- 
lations are performed at half filling (n) = 1. In the har- 
monic case, we have fi = 0, but the anharmonic problem 
must have the chemical potential adjusted as a function 
of temperature to yield the right filling. 

When the system is in an ordered charge-density-wave 
state, there are two probability distributions for the 
phonon coordinate — one for the A sublattice wa{x) and 
one for the B sublattice wb(x). The order parameter for 
the charge-density- wave phase is defined to be 



A(T) = g / dx[w A (x) - w B (x)]x, 



(3) 



which measures the average difference in the phonon co- 
ordinate between the A and B sublattices multiplied by 
the deformation potential. 
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FIG. 1. Charge-density- wave order parameter as a function 
of temperature for three cases of anharmonicity (all chosen to 
be at the maximal T c ) (i) g = 2, a an = 0.16, T c = 0.1277 
(solid line); (ii) g = 1.5, a a „ = 0.03, T c = 0.1326 (dashed 
line); and (iii) g = 1.19, a a „ = 0, T c = 0.1340 (dotted line). 
The figure is renormalized to plot A(T) / A max and T/T c . In- 
set is a plot of the average phonon coordinate on the A and 
B sublattices. Note that all the curves are symmetric, but 
their average value moves away from zero for the anharmonic 
cases. 

In Figure 1 we show a plot of this charge-density-wave 
order parameter as a function of temperature. The plots 
include three cases, all at the maximum value of T c for 
a given value of g and arbitrary a an . The parameters 
chosen are g = 2, a an — 0.16 which has T c = 0.1277; 
g = 1.5, a an = 0.03 which has T c = 0.1326; and g = 1.19, 
a an = which has T c = 0.1340. The curves are normal- 
ized by the maximal gap value calculated at T w T c /8 
which is very close to the T = value, and by T c . Note 
that rather than approaching 1, there is a small decrease 
in A(T) / A max at the lowest temperatures. This is an ac- 
curacy issue with our calculations (which are performed 
with a fixed number of 5000 Matsubara frequencies for 
the energy cutoff). We estimate that our error in the 
T — gap is no larger than about 1%. These curves have 
the correct generic behavior that we expect for mean- 
field systems: the order parameter increases like \JT C — T 
away from the transition temperature and then rapidly 
saturates. We see no significant difference in the shape 
of these curves for different values of the anharmonicity. 
Inset into this figure is a plot of (xa) upper curves and 
(xb) lower curves. Note how the anharmonic systems are 
both shifted upwards and have narrower spreads in the 
average phonon coordinate, but that the curves are sym- 
metric between the A and B sublattices, so the particle- 
hole asymmetry only affects the midline of the phonon 
coordinate, not its distortion in the ordered phase! 

In summary, we can accurately estimate the value of 
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the T = gap by performing calculations at T — T c /10 
and the shape of the order parameter (as a function of 
T) is not too strongly dependent on the strength of the 
anharmonicity. In addition, we find that the main ef- 
fect of the anharmonic interaction is to change the aver- 
age of the phonon coordinate values on each sublattice 
{(xa) + (xs))/2 and shrink the magnitude of the dis- 
tortion. This does not mean that the gap is reduced 
by as much, though, because systems that share the 
same approximate value of T c will have different values 
of both g and a an . The reduction in the distortion of the 
phonon coordinate can be compensated by a correspond- 
ingly larger value of g. 



III. RESULTS AND DISCUSSION 

We now present results for the charge-density- wave gap 
as a function of the anharmonicity. The most reasonable 
way to present these results is to plot the gap versus a 
measure of the anharmonic potential energy in equilib- 
rium. We can determine what the equilibrium phonon 
coordinate x — x* is for the atomic problem [determined 
by Eq. ([!]) with — and (n) = 1] and then plot results 
versus a Qn a;* 4 . These results are summarized in Fig. 2. 
Note that this measure of the anharmonic potential en- 
ergy is double valued when g is large, since it approaches 
in both the small and large a an limits. The general 
shapes of these curves are quite similar, but the scales 
change with the coupling strength and the stronger cou- 
pled cases show more curvature (and eventually a double- 
valuedness) . 
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FIG. 2. Charge-density- wave gap calculated at T c /10 plot- 
ted versus the anharmonic "energy" a an x* 4 determined from 
the average value of the phonon coordinate for the atomic 
problem. Note how the different curves (g = 0.725, 0.95, 
1.19, 1.5, and 2.0) have the same general shape, which be- 
come double valued for the strongest coupling cases shown 
here. 



We can examine the static Holstein -model in both a 
strong-couplingcJ and a weak-couplingi3 limit. In the 
strong coupling limit, the transition temperature ap- 
proaches zero as l/(2g[xo — X2]), while the zero temper- 
ature gap becomes large as g(xa — x-i)j1 and hence the 
ratio can diverge [2A/(/cbT c ) = 2g 2 (xo — X2) 2 ]. Here Xi 
satisfies 



g(i - 1) + nxi + 3f3 an x 2 + 4cv Q 



0. 



(4) 



This ratio becomes infinite as the coupling strength in- 
creases, hence there is no limit to the magnitude of the 
gap ratio. Note that this result is also true for the su- 
perconducting order, when one includes the effect of a 
nonconstant density of states, since the strong-coupled 
T c will approach zero there as well. This result comes en- 
tirely from the fact that the phonon-coordinate distortion 
grows linearly with g in strong coupling, but the transi- 
tion temperature decreases as 1/g 2 due to the strong- 
coupling physics. 

We can also investigate the possibility of universal be- 
havior in strong coupling. Since the system consists of 
empty sites and preformed pairs, the self energy has a 
low-energy pole and takes the form S(w) = 01/ (u> + iS) + 
0(u> ) for small \ui\ with a > 0. Plugging this form 
into the self-consistent equations for the Green's function 
yields a = — 1/2 + g 2 \xo\\x2\ (which is larger than zero 
for g large enough). Defining the wavefunction renor- 
malization parameter by a scaling along the imaginary 
axisliil 
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with iujj = inT(2j + 1), yields 



Z = 1 



13g 2 \x \\x 2 \ 
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(5) 



(G) 



Since Z is a function of xq and X2 and 2A/(fcsT c ) is a 
function of X0 — X2, we don't expect universal behavior at 
extremely strong coupling for the anharmonic case. But, 
in the harmonic case we have |a;o| = \%2\, so we expect 
deviations from universality only when |xo|/|x2| deviates 
far from unity. 

The weak-coupling limit is much more co mp l ic ated. 
There are many approaches that can be takenBliU, but 
none produce good agreement with the transition tem- 
perature over a wide range of coupling strengths (but the 
zero-temperature gap is approximated well) . Here we will 
concentrate on just two different strategies (for the ha*-i 
monic case only): (i) the renormalized phonon methodtj 
where a certain class of vertex correction terms can be 
neglected from the analysis, but one needs to work with 
renormalized phonons (renormalized by the electroa-holc 
bubble diagrams) and (ii) a similar approximation^ that 
employs the identical set of diagrams, but does not renor- 
malize the phonons. 
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FIG. 3. Ratio of twice the charge-density-wave gap (calcu- 
lated at T c /10) to the charge-density- wave transition temper- 
ature T c plotted versus the wavefunction renormalization pa- 
rameter. The harmonic limit is connected by a thin solid line, 
while the renormalized (dotted) and unrenormalized (dashed) 
weak coupling approximations are also included. Note how 
the anharmonic systems generically lie close to the harmonic 
results for a wide range of parameters. Deviations occur 
for large g and for intermediate values of the net coupling 
strength, but generally are no larger than about 20% (except 
at very strong coupling where the gap ratio is enhanced). 

In Fig. 3 we show the results for the gap ratio in the 
harmonic case, five different anharmonic cases, and the 
two approximation schemes described above. In order to 
check for universal behavior, we plot the results as a func- 
tion of the wavefunction renormalization parametertil Z 
[which is evaluated as in Eq. (^)]. One can see some 
striking behavior in this plot. First, we find that the gap 
ratio can become as large as one would like as the cou- 
pling strength gets bigger and bigger. Furthermore, the 
gap ratio rises very rapidly above the weak-coupling limit 
of 3.52, so even Z factors of 1.25 (which would correspond 
to very weak electron-phonon coupling) would have a gap 
ratio larger than 6. We find that the calculated results lie 
in between the two different weak coupling schemes in- 
dicating that phonon renormalization is important, but 
the summation of the bubble diagrams renormalizes the 
phonons too strongly in the general case. Finally, we 
notice that except for the case of very strong coupling 
(g = 2), the gap ratio is modified by at most about 20% 
from the harmonic case. This finding is most surpris- 
ing, but can be explained in relatively simple terms. As 
the anharmonic potential increases, the system becomes 
weaker and weaker coupled, because the anharmonic po- 
tential prevents the phonon coordinate from deviating 
far from the origin. In the limit of very weak coupling, 
the transition is dominated by a nesting instability and 
can be described within a BCS-like format. Hence all 
results must agree in this limit (if one properly identifies 
the net strength of the attractive interaction) . Further- 



more, when the system has small anharmonicity, the an- 
harmonic effects can be treated perturbatively, and the 
system remains close to the harmonic limit there as well. 
Since the curves are pinned to be close to the harmonic 
limit for small anharmonicity and for large anharmonic- 
ity, we find that they generically do not stray far from the 
harmonic curve in the intermediate regime. As the cou- 
pling strength increases so this intermediate regime be- 
comes larger, the deviations can also become significant, 
as we see for the g = 2.0 case. But we find it surprising 
that there is such a wide range of parameter space where 
the results for the anharmonic system remain so close to 
the harmonic system! 

IV. CONCLUSIONS 

Contrary to the simple arguments about anharmonic- 
ity, we find that generically anharmonic phonons are de- 
scribed well by an equivalent harmonic limit, even in 
the ordered phase, where one would expect the effects 
of the anharmonicity to be felt more strongly. Hence, 
we find that the analysis given for the large gap ratio 
in the harmonic electron-phonon problem remains essen- 
tially unchanged as anharmonicity is introduced except 
in the limit of extremely strong electron-phonon coupling 
(which may be so unphysically large that is is not at- 
tained in any real material). The reason why this holds 
is essentially a continuity argument: for small and large 
anharmonicity, the system must be close to the harmonic 
results — hence it remains close for intermediate values as 
well. |— - 

These results imply, once againEJ, that the quasihar- 
monic approximation should be quite accurate in real ma- 
terials, since the anharmonic system can be mapped onto 
an effective harmonic system that shares both the single- 
particle and two-particle properties. Deviations are ex- 
pected only in the strong coupling and strong anharmonic 
limit. The deviation of the gap ratio from the BGS value 
can be understood from the harmonic analysiaij. We 
don't expect these results to be changed much for small 
phonon frequencies, because we know that T c follows a 
universal form for low frequency, and we expect the gap 
value will not change much if the phonon frequency is 
small. 
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